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To gain deeper insight into the dynamics of complex quantum systems we need a quantum leap 
in computer simulations. We can not translate quantum behaviour arising with superposition states 
or entanglement efficiently into the classical language of conventional computers. The final solution 
to this problem is a universal quantum computer If], suggested in f982 and envisioned to become 
functional within the next decade(s); a shortcut was proposed via simulating the quantum behaviour 
of interest in a different quantum system, where all parameters and interactions can be controlled 
and the outcome detected sufficiently well. 

Here we study the feasibility of a quantum simulator based on trapped ions We experimentally 
simulate the adiabatic evolution of the smallest non-trivial spin system from the paramagnetic into 
the (anti-)ferromagnetic order with a quantum magnetisation for two spins of 98%, controlling 
and manipulating all relevant parameters of the Hamiltonian independently via electromagnetic 
fields. We prove that the observed transition is not driven by thermal fluctuations, but of quantum 
mechanical origin, the source of quantum fluctuations in quantum phase transitions "3]. We observe 
a final superposition state of the two degenerate spin configurations for the ferromagnetic (|Tt) + UJ-)) 
and the anti-ferromagnetic (|T!) + |!T)) order, respectively. These correspond to deterministically 
entangled states achieved with a fidelity up to 88%. 

Our work demonstrates a building block for simulating quantum spin-Hamiltonians with trapped 
ions. The method has potential for scaling to a higher number of coupled spins [2|. 

PACS numbers: 



I. INTRODUCTION 

It is not possible to efficiently describe the time evo- 
lution of quantum systems on a classical device, like a 
conventional computer, since their memory requirements 
grow exponentially with their size. For example, a clas- 
sical memory needs to hold 2 50 numbers to store arbi- 
trary quantum states of 50 spin-1/2 particles. To be 
able to calculate its evolution demands to derive a ma- 
trix of (2 50 ) 2 = 2 100 elements, already exceeding by far 
the capacity of state of the art computers. Each dou- 
bling of computational power permits only one additional 
spin-1/2 particle to be simulated. To allow for deeper in- 
sight into quantum dynamics we need a "quantum leap" 
in simulation efficiency. 

As proposed by Richard Feynman [l| , a universal quan- 
tum computer would accomplish this step. A huge va- 
riety of possible systems are under investigation, a very 
promising one being trapped ions [1, [f| acting as quan- 
tum bits (qubits). After addressing the established crite- 
ria summarized by DiVincenzo @ on up to 8 ions 0, [1] 
with operational fidelities exceeding 99% 0, H, Q , there 
seems to be no fundamental reason why such a device 
would not be realisable. 

An analogue quantum computer, much closer to the 
original proposal by Feynman, might allow for a short- 
cut towards quantum simulations. We want to simulate 
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a system by a different one being described by a Hamilto- 
nian that contains all important features of the original 
system. The simulator needs to be controlled, manipu- 
lated and measured in a sufficiently precise manner and 
has to be rich enough to address interesting questions 
about the original system. For large coupled spin systems 
optical lattices might be advantageous [l(| , while smaller 
spin systems and degenerate quantum gases might be 
simulated by trapped ions 0, QjJ . Instead of implement- 
ing a Hamiltonian with a universal set of gates, direct 
simulation of the Hamiltonian typically consists of one 
(adiabatic) evolution of the initial state into the corre- 
sponding final state of interest. 

Here, in a proof-of-principlc experiment, we simulate 
the adiabatic transition from a quantum para- to a quan- 
tum (anti-)ferromagnet and illustrate the advantages of 
the adiabatic quantum simulation (see FIG. [J). We 
demonstrate the individual access, via rf- and laser fields, 
to all relevant parameters in the underlying Hamiltonian, 
representing one out of a large spectrum of quantum spin- 
Hamiltonians. 



II. ADIABATIC QUANTUM SIMULATION 

The adiabatic quantum simulation of generic spin- 
Hamiltonians proposed by Porras and Cirac [2| can be 
illustrated considering a string of charged spin-f/2 par- 
ticles confined in a common harmonic potential. Two 
electronic states of each ion simulate the two-level sys- 
tem of a spin-f/2 magnetic moment, ||) and |J,). Note 





FIG. 1: Phase transition of a quantum magnet: Each 
ion can simulate a magnetic spin, analogue to an elemen- 
tary magnet. We initialise the spins in the paramagnetic 
state |— »— > . . . — >), the ground state of the Hamiltonian 
Hb = B x (af + erf + • • • + This is equivalent to align- 

ing the spins parallel to the simulated magnetic field. Adding 
an effective spin-spin interaction J(t) (at constant B x ) and 
increasing it adiabatically to | J m ax | 3> B x , we expect the sys- 
tem to undergo a quantum phase transition into a ferromag- 
net, the new ground state of the system (J is symbolized as lit- 
tle chains, trying to keep neighbouring spins aligned) . Ideally, 
the two possible ferromagnetic orders |TT . . . |) and 111 . . . |) 
are degenerate ground states. The spin system should evolve 
into the superposition state |TT ••• T) + 1 II ••• l)> a maximally 
entangled "Schrodinger Cat" state/magnet. 



FIG. 2: Simulating the quantum magnet: a) Two perpendicu- 
lar polarized laser beams of frequency wi and LU2 induce a state 
dependent optical dipole force F± = — 3/2 i*\ along the trap 
axis a by the AC-Stark shift (here, only F-\ is depicted), b) 
For a standing wave u)i = W2, the force conditionally changes 
the distance between neighbouring spins simulating a spin- 
spin interaction [2J (F^ (F\) symbolised by the arrow to the 
right (left)). Only if all spins are aligned (top), the total 
Coulomb energy of the spin system is not increased, defining 
ferromagnetic order, the quantum magnet, to be the ground 
state. For u>i 7^ L02, the sinusoidal force pattern can be seen 
as a wave moving along the trap axis a pushing or pulling the 
ions repeatedly at a frequency wi — ui2- We chose U\ — U02 
close to the resonance frequency of the ions oscillating out of 
phase (stretch mode). The energies of different spin states 
now depend on the coupling of the spin state to the stretch 
mode. Energy can be coupled efficiently into the state with 
different spin orientations (e.g. bottom), defining the not af- 
fected upper case as a ground state [l4|]. The interpretation in 
terms of an effective spin-spin interaction is further described 
in methods. 



that the inter-ion distance of several um renders any di- 
rect spin-spin coupling negligible. The quantum Ising 
Hamiltonian, 

Rising = H B + Hj = B x of + J v a i a h W 

i i<j 

consists of two terms. The first denotes the interaction of 
each individual spin, represented by the Pauli operator 
<j\ (k can be x, y, or z), with a uniform magnetic field 
B x pointing into direction x. The second term stands for 
the spin-spin interaction which tries to align the spins 
(erf) parallel or anti-parallel along the z-axis dependent 
on the sign of the interaction amplitude Jy . To simulate 
the first one they couple the eigenstates of erf, ||) and 
||), via an electromagnetic field. The latter is simulated 
by a state-dependent forcing [12, G3 > further explained 
with the help of FIG. H 

To understand the experiment discussed below, we can 
consider two extreme scenarios and interactions between 
nearest neighbours only, Ji.i+i = J. For the case of J = 
and B x > 0, the ground state of the spin-system has all 
spins aligned with B x along the x-axis. This corresponds 
to the paramagnetically ordered state |— ►— * . . . — ►), the 
eigenstate of the Hamiltonian B x ^2 i af with the lowest 
energy. 

For the opposite case of B x — and J < (J > 0), 
the system has an infinite number of degenerate ground 



states, defined by any superposition of the lowest 
energy eigenstates of J2i<j J y°f cr l ■> namely |TT---t) 
and HI • • • I) which represents ferromagnetic order (or 
ITIT • • • Tl) an d IITI • • • IT), the anti-ferromagnetic or- 
der, respectively). Initialising the spin system in an 
eigenstate in the cr^-basis, starting with J(t = 0) = 
and B x > and adiabatically increasing |J(t)| to 
I Jmax I ^ B x should evolve the system from the param- 
agnetic arbitrarily close into the (anti-)ferromagnetic or- 
der, as depicted in FIG. [TJ A quantum phase transition 
is supposed to occur at B x = \J\ in the thermodynamic 
limit of an infinite amount of spins [H, [l5[ • 



III. EXPERIMENTAL IMPLEMENTATION 

We implemented the experiment as described in the 
following. We confine two 25 Mg + ions in a linear Paul 
trap [16[ and laser-cool them to the Coulomb-crystalline 
phase where the ions align along the trap axis a. The 
motion of the ions along a can be described in the basis of 
normal modes: the oscillation-in-phase-mode (com) and 
the oscillation-out-of-phase-mode (stretch). The related 
oscillation frequencies amount to w com = 2ir x 2.1 MHz 
and w s tretch = 27r x 3.7 MHz, respectively. 

In our implementation we define the hyperfine ground 
states ||) = \F = 3, m f = 3) and |f) = \F = 2,m f = 2) 
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in the 2S1/2 levels separated by ui = 2tt x 1.7 GHz. An 
external magnetic field B of 5.5 G orients the magneti- 
sation axis for the projection hrrif of each ion's angular 
momentum F. In this field adjacent Zceman sublevels of 
the F = 3 and F = 2 manifolds split by 2.7 MHz per 
level. 

We coherently couple the states ||) and the ||) with a 
resonant radio-frequency field at loq to implement single 
spin rotations 0, 113] , 

R{e, <f) = cos(6/2) / - i sin(9/2) cos(0) a x 

-ism(&/2)sm(4>)a y , ' 2 ' 

where / is the identity operator, a x and a v denote the 
Pauli spin matrices acting on |1) and ||), 9/2 = B x t is 
proportional to the duration t of the rotation and <f> is the 
phase of the rf-oscillation, defining the axis of rotation in 
the x-y-plane of the Bloch sphere. 

We provide the effective spin-spin interaction by a state 
dependent optical dipole force [ID, [TJ, [l8j . The relative 
amplitudes F± = — 3/2 are due to AC-Stark shifts in- 
duced by two laser beams at wavelength A of 280 nm, 
depicted in FIG. 2a, perpendicular in direction and polar- 
isation with their effective wave-vector difference point- 
ing along the trap axis a. They are detuned 80 GHz 
blue of the 2 P 3 / 2 excited state, with intensities allow- 
ing J/H above 2tt x 22.1 kHz. We use a walking wave 
force-pattern by detuning the two laser wavelengths by 
2?r x 3.45 MHz = w str ctch + 6 with S = -2tt x 250 kHz. 
This choice avoids several technical problems of the orig- 
inal proposal @ (see methods), while at the same time, 
resonantly enhancing the effective spin-spin interaction 
by a factor of [kWetch/^1 = 14.8 compared to the stand- 
ing wave case [191 ] . 

After laser cooling we initialise the quantum simula- 
tor by optical pumping Q to the state |J.)|J.) \n = 0). 
We rotate both spins in a superposition state via a 
R(tv/2, — 7r/2)-pulse (see Eq. 2) on the rf-transition to ini- 
tialise the state = |— = 0). Note that this 
paramagnetic state |— ►>]—►> = (|T) + |I))(|T) + ID) = 
ITT) + ITI) + I IT) + III) has a 25% probability to be pro- 
jected into either |TT) or |J.J.) (normalisation factors are 
suppressed throughout). 

After the adiabatic evolution described below, we 
project the final spin state into our a z -measurement 
basis by a laser beam tuned resonantly to the 
|T) <-> 2 P 3 / 2 \F = 4, to/ = 4) cycling transition Q. An 
ion in state ||) fluoresces brightly, leading to the detec- 
tion of on average 40 photons during a 160 us detection 
period with our photo multiplier tube. In contrast, an ion 
in state |T) remains close to dark (on average 6 photons). 
We repeat each experiment for the same set of parameters 
10 4 times and derive the probabilities Pn,P^ and P^ 
for the final state being projected into state |||) , |TT) 
and ||T) or |T|), respectively (and further described in 
methods) . 

We simulate the effective magnetic field by continu- 
ously applying a radio-frequency field with phase = 
and an amplitude such that it corresponds to a single 




FIG. 3: Quantum magnetisation of the spin system: We 
initialise the spins in the paramagnetic state |— >) |— ») = 
(IT) + IDXIT) + ID) = ITT) + ITI) + III) + III), the ground 
state of the Hamiltonian Hb = B x (af +erf). A measurement 
of this superposition state would already project into each 
|tt) and HI) with a probability of 0.25. After applying B x 
we adiabatically increase the effective spin-spin interaction 
J(t — 0) = to J(T). State sensitive fluorescence detec- 
tion allows to distinguish the final states |TT), III), III) or 
lit). Averaging over 10 4 experiments provides us with its 
probability distribution Pn (two ions fluoresce), Pj-j (no ions 
fluoresces), and Pjj or Plj (one ion fluoresces). We repeat 
the measurement for increasing ratios J(T)/B X . The exper- 
imental results for the ferromagnetic contributions Pff and 
Pi I are depicted as squares, the solid lines representing the 
theoretical prediction. For J(T)/B X <C 1, the paramagnetic 
order is preserved. For J(T)/B X 3> 1, the spins undergo a 
transition into the ferromagnetic order, the ground state of 
the Hamiltonian Hj = J max o"iO"2, with a related quantum 
magnetisation M = Pj.j. + P TT of > (98 ± 2)%. Note that 
we invert the ordinate of the lower frame to emphasise the 
unbroken symmetry of the evolution. 



qubit rotation i?(9,0) with full rotation period G = 2ir 
in 118 us and deduce B x = 2ir x 4.24 kHz. Precise control 
of the phase (f> of the rf-oscillator relative to the initiali- 
sation pulse allows to align B x parallel to the spins along 
the s-axis in the equatorial plane of the Bloch sphere, en- 
suring that is an eigenstate of this effective magnetic 
field. 

At the same time, we switch on the effective spin-spin 
interaction J(i) (t £ [0;T]) and increase its amplitude 
adiabatically up to J(T). At time T, we switch off the 
interactions and analyse the final state of the two spins 
via the state sensitive detection described above. In a 
sequence of experiments at constant B x we increase T 
and therefore J(t)/B x . After 50 steps of 2.5 ps each 
we reach the maximal amplitude J(t = 125 ]xs)/B x = 
Jmax/B x = 5.2 (see methods) and achieve a quantum- 
magnetisation M. = Pi i +-Ptt, the probability of being 
in a state with ferromagnetic order, of M. = (98 ± 2)%. 
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FIG. 4: Entanglement of the quantum magnet: Measurement 
of the parity V = Pn + Pn — (Pit + P\l) °f the nna l ferromag- 
netic state after the simulation reached J m a*/B x — 5.2. As we 
vary the phase of a subsequent analysis pulse, the parity of 
the two spins oscillates as C cos(2<^>). Together with the final 
state populations Pjj and P-fj depicted in FIG. [3] we can de- 
duce a lower bound for the fidelity T = 1/2(P U + P TT ) + C/2 
of (88 ± 3)% for the final superposition state |ff) + Iff), a 
maximally entangled state, highlighting the quantum nature 
of this transition. We find qualitatively comparable results 
for the antiferromagnetic case |ff) + |ff). Each data point 
averages 10 4 experiments. 



IV. ENTANGLEMENT 

In our experiment we can detect both ferromagnetic 
contributions, Pjj and Pf^, separately. Any imperfec- 
tion in the simulation acting as a bias field B z along the 
z-axis, would energetically prefer one of the ferromag- 
netic states over the other and therefore unbalance their 
contribution to the final state. We carefully cancel all 
bias fields (see methods) to balance the populations Pf j 
and Pn, as can be seen in FIG. [3] The results are in 
good agreement with theoretical predictions for our ex- 
periment, shown as solid lines. We expect the final state 
to be a coherent superposition of the two ferromagnetic 
states |ff) + |ff), close to a maximally entangled Bell 
state. To quantify the experimentally reached coherence 
we measure the parity [21| V = Pj.x + Pn — (-P.it + -Pn) 
after applying an additional R(n/2, </>)-pulse to both ions 
after J max is reached, with a variable rf-phase (f> relative 
to the rf-field simulating B x . The measured data shown 
in FIG. 0] have a component that oscillates as C cos(2</>), 
where |C|/2 characterises the coherences between the |f f ) 
and |ff) components in the state produced. Deducing a 
contrast C of (78±2)% from the best-fit we derive a lower 
bound of the fidelity [p T = 1/2(P U + P TT ) + C/2 of 
(88 ±3)%. 

We also simulate the adiabatic evolution of a sys- 
tem not initialised in the ground state of the initial 



Hamiltonian. In particular we prepare the paramagnetic 
eigenstate |<— <— ) = (||) - |T))(|t) - ID); witu the spins 
aligned anti-parallel with respect to the simulated mag- 
netic field via a R(tt/2,tt/2) rf-initialisation pulse. The 
adiabatic evolution should preserve the spin system in 
its excited state leading now into the anti-ferromagnetic 
order |ff) + |ff). After evolution to J — J max we 
find P iT + P Ti > (95 ± 2)%. To investigate the coher- 
ence between the |ff) and the |ff) components, we ro- 
tate the state via an additional R(tt/2, 0)-pulse which 
would ideally take |f|) + |ff) — > |ff) + |f|), before 
we continue to measure the parity as explained above. 
We deduce a lower bound for the fidelity of the anti- 
ferromagnetic entangled state T — | (^fi n ai | if + f |) | = 
l/2(P n +P iT ) + C/2of (80 ±4)%. 

An equally valid viewpoint of this experiment inter- 
prets |<— <— ) as the ground state of the Hamiltonian 
—Poising. Because the sign of all spin-spin interactions 
is also reversed in —Poising it is equivalent to a change of 
sign in the spin-spin interaction J. 

The entanglement of the final states additionally con- 
firms that the transition from paramagnetic to (anti-) 
ferromagnetic order is not caused by thermal fluctua- 
tions driving thermal phase transitions. The evolution 
is coherent and quantum mechanical, the coherent coun- 
terpart to the so-called quantum fluctuations [!, driv- 
ing quantum phase transitions in the thermodynamic 
limit. In this picture tunnelling processes [l5| induced 
by B x coherently couple the degenerate (in the rotating 
frame) states |f) and |f) with an amplitude proportional 
to (B x /\ J\). In a simplified picture for N spins the ampli- 
tude for the tunnelling process between ^ m] = If f • ■ ■ f) 
and 'Fjvi = Iff • • ■ I) is proportional to (B X /\J\) N , since 
all N spins must be flipped. In the thermodynamic limit 
(N — > oo) the system is predicted to undergo a quan- 
tum phase transition at |J| = B x . At values J > B x 
the tunnelling between 'J'ocj and 'Fool is completely sup- 
pressed. In our case of a finite system an d ^21 re- 
main coupled and the sharp quantum phase transition 
is smoothed into a gradual change from paramagnetic to 
(anti-)ferromagnetic order. 



V. CONCLUSION AND OUTLOOK 

We demonstrated the feasibility of simple quantum 
simulations in an ion trap by implementing the Hamil- 
tonian of a quantum magnet undergoing a robust transi- 
tion from a paramagnetic to an entangled ferromagnetic 
or anti-ferromagnetic order. While our system is cur- 
rently too small to solve classically intractable problems, 
it uses an approach that is complementary to a universal 
quantum computer in a way that can become advanta- 
geous as the approach is scaled to larger systems. Since 
our scheme only requires inducing the same overall spin- 
dependent optical force on all the ions , it does not rely 
on the use of sequences of quantum gates, thus its scaling 
to a higher number of ions can be simpler. Furthermore 
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the desired outcome might not be affected by decoherence 
as drastically as typical quantum algorithms, because a 
continuous loss of quantum fidelity might not spoil com- 
pletely the outcome of the experiment (for example the 
(anti-)ferromagnetic ordering transitions are hardly af- 
fected by phase decoherence), while universal quantum 
computation will almost certainly require involved sub- 
agorithms for error correction [J]. Decoherence in the 
simulator might even mimic the influence of the natural 
environment [22j | of the studied system, if we judiciously 
construct our simulation (for example the decoherence 
we mainly observe in our demonstration implements a 
dephasing environment). 

Despite technical challenges, we expect that this work 
is the start to extensive experimental research of com- 
plex many-body phases with trapped ion systems. Lin- 
ear trapping setups may be used for the quantum simula- 
tion of quantum dynamics beyond the ground state where 
chains of 30 spins would already allow to outperform cur- 
rent simulations with classical computers. We may also 
adapt our scheme to new ion trapping technologies [l3|. 
For example, a modest scaling to systems of 20 x 20 spins 
in 2D would yield insight into open problems in solid state 
physics, e.g. related to spin-frustration. This could pave 
the way to address a broad range of fundamental issues 
in condensed matter physics which are intractable with 
exact numerical methods, like, for example, spin liquids 
in triangular lattices, suspected to be closely related to 
phases of high-T c superconductors [H[ . 



VI. METHODS 



the axial trapping frequency to better than 100 Hz. In 
contrast to phase gate implementations [l8[ we have to 
detune the two laser beams far enough to keep the mo- 
tional excitation and the related errors due to residual 
spin-motion coupling [2j insignificant. Adjusting the de- 
tuning to 5 = — (wgtretch - (^l — ^2)) = — 27T x 250 kHz red 
of the stretch mode frequency and terminating the evo- 
lution after the system returned back into the motional 
ground state [l8| ideally completely cancels the simula- 
tion errors discussed in @. 

State sensitive detection: For two spins, the integrated 
fluorescence signal does not allow to distinguish between 
two (|tl) and |||)) of the four possible spin configura- 
tions. In addition, the amount of detected photons for 
each of the three distinguishable configurations fluctuates 
from experiment to experiment according to Poissonian 
statistics and therefore can only be determined with lim- 
ited accuracy. For the data reported, we repeated each 
experiment 10 4 times and fitted the resulting photon- 
number distribution to the weighted sum of three refer- 
ence distributions to derive Pn, P^ and Pfj + Pit- 

Adiabatic evolution: We achieve the best fidelities for 
the reported transitions at a duration of the simulation 
of T = 125 ps at a B x — 2tt x 4.24 kHz. Even though we 
are not strictly in the adiabatic limit, the robustness of 
the transition allows to minimise decoherence effects re- 
ducing the duration of the simulation. In addition, tech- 
nical reasons led to the evolution of J(t) linear in time 
to up to J(t = 50 ps) = 5 x 10~ 4 J ma x, continued by 
J(t) « (e Qt - (3f best fitted by a = 0.026 and = 4. 
Up to now we did not improve the fidelities by evolving 
or terminating J(t) or B x in a more adiabatic way. 



State dependent optical dipole force: An effective 
(Ising) spin-spin interaction was proposed to be im- 
plemented via magnetic field gradients [13]. Por- 
ras and Cirac su gge sted to use state dependent opti- 
cal dipole forces [ill [l8j displacing the spin state \s) 
(s either f or 1) in phase space by an amount that de- 
pends on \s). The area swept in phase space changes 
the state to e^ s ' \s). The phase <p{s) can be broken 
down into single spin terms proportional to of and ap- 
parent spin-spin interactions proportional to crfaj and 
thus gives rise to the desired simulation of spin-spin in- 
teractions 25]. It can also lead to single spin phases 
that simulate the unwanted contribution of a common 
bias fields B z a z in the Hamiltonian that will unbalance 
the probabilities Pyi and Pff. To achieve a balanced 
probability distribution as depicted in FIG. 3, we have 
to carefully compensate these single spin phases. To this 
end we (1) compensate the residual AC-Stark shifts of 
the individual laser beams by carefully choosing direc- 
tion and polarisation of the beams [l(| and (2) compen- 
sate for the imbalance caused by single spin phases via a 
detuning of the order of several kHz of the rf-transition 
relative to ujq. (3) The ions have to be separated by an 
integer multiple of the effective wavelength A e ff = A/v2, 
in our implementation 18 x A e ff, requiring the control of 
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